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Abstract 

We apply the S-matrix based finite temperature formalism to non-relativistic Bose and Fermi 
gases in 1 + 1 and 2 + 1 dimensions. For the 2 + 1 dimensional case in the constant scattering 
length approximation, the free energy is given in terms of Roger's dilogarithm in a way analagous 
to the thermodynamic Bethe ansatz for the relativistic 1 + 1 dimensional case. The Id fermionic 
case with a quasi-periodic 2-body potential could provide a physical framework for understanding 
the Riemann hypothesis. 
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I. INTRODUCTION 



Quantum gases at finite temperature and density have very wide applications in nearly 
all areas of physics, ranging from black body radiation, Fermi liquids, Bose-Einstein conden- 
sation, and equations of state in cosmology. If one knows the zero temperature dynamics, 
i.e. the complete spectrum of the hamiltonian, then the quantum statistical mechanics just 
requires an additional statistical summation Z = Tr(e _/3H ), so that in principle zero tem- 
perature dynamics and quantum statistical sums are decoupled. This decoupling is also 
in principle clear from an intuitive picture of a gas as a finite density of particles that are 
subject to scattering. In practice, the complete non-perturbative spectrum of H is unknown 
so one must resort to perturbative methods, such as the Matsubara approach, that typically 
entangle the zero temperature dynamics from the quantum statistical mechanics. 

In |l| an alternative approach to finite temperature quantum field theory was developed 
that achieves this decoupling of zero temperature dynamics and quantum statistical sum- 
mation. The dynamical variables are the occupation numbers, or filling fractions, of the 
gas, and the zero temperature data of the underlying theory are just the S-matrix scattering 
amplitudes. Our construction was modeled after Yang and Yang's thermodynamic Bethe 
ansatz (TBA)Q). The TBA is very specific to integrable theories in one spacial dimension. 
Nevertheless, our derivation [lj is carried out in any spacial dimension d and doesn't assume 
integrability; therefore it is in general an approximate method, though it can be systemat- 
ically improved. An important ingredient of our construction is the work of Dashen, Ma, 
and Bernstein [3j which explains how to formulate quantum statistical mechanics in terms of 
the S-matrix. Ideas of Lee and Yang were also instrumental^], even though the latter work 
is not based on the S-matrix. 

Our previous work[l| was mainly devoted to developing the formalism in a general, model 
independent way for both non-relativitistic and relativistic theories. Though the derivation 
was somewhat involved, the final result is straightforward to implement and is summarized 
by the two formulas ( 1231 \2M . In this paper we apply the method to non-relativistic Bose 
and Fermi gases with special attention paid to the low dimensional cases of d = 1, 2. The 
main approximation we make is to consider only the contributions to the free energy that 
come from 2-particle to 2-particle scattering. On physical grounds these are expected to be 
the most important if the gas is not too dense. Even in this approximation the problem 
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is non-trivial since one needs to resum all the 2-body interactions self-consistently. This is 
accomplished by an integral equation for the filling fractions that is analagous to the TBA 
equations. 



In the next section we first review the main results in |l|. In section III we describe the 
kind of second-quantized hamiltonians that are the subject of this paper. In section IV we 
consider bosons with a 5 function two-body potential. In section V we turn to fermionic 
gases in the same approximation and derive formulas that can be used to study the effect 
of interactions on the Fermi energy. Since the lowest order contribution corresponds to a 
constant scattering length approximation, the results we obtain in sections IV and V are 
already known, and are included here mainly for illustrative purposes. 

In section IV we study the 2d case which has some remarkable features. The integral 
equation which determines the filling fraction becomes algebraic. The formulas for the free 
energy are essentially identical to those for the conformally invariant limit of relativistic 
theories in one lower dimension d = 1 [5| and are given in terms of Roger's dilogarithm. In 
making this analogy we define a "central charge" c as the coefficient in the free energy, which 
in the relativistic case in Id is the Virasoro central charge. As in the Id relativistic case, for 
certain special values of the coupling, c can be a rational number, and some examples based 
on the golden mean are presented. We also extend the formalism to many species of mixed 
bosonic and fermionic statistics. 

As will be clear in the sequel, the polylogarithmic and Riemann ( functions play a cen- 
tral role in this work. This provided us with a new opportunity to understand Riemann's 
hypothesis [fj. Riemann's ( function for Re{v) > 1 is defined as 



It has a simple pole at v — 1; as we will describe this pole is a manifestation of the impos- 
sibility of Bose-Einstein condensation in 2d. The function can be analytically continued to 
the whole complex v plane. It has trivial zeros at v = —2n where n is a non-zero positive 
integer. The function has a duality relation that relates ((u) to ((1 — u) (eq. (11091) below) 
so that its non-trivial zeros are in the critical strip < Re{p) < 1 and symmetric around 
Re(u) = 1/2. The Riemann hypothesis (RH) is the statement that the only non-trivial zeros 
of ((u) are at Re(u) = 1/2. This hypothesis is very important in number theory since, as 
shown by Riemann, the distribution of prime numbers is intimately related to the location 
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of the Riemann zeros. It has already been proven by Hadamard and de la Vallee-Poussin 
that there are no zeros with Re{v) = 1, which is equivalent to proving the Prime Number 
Theorem 6| . 

There have been a number of approaches to understanding the Riemann hypothesis based 
on physics 26|]. Here we mention a few of them mainly to contrast with the approach pre- 
sented here. Some approaches are inspired by the Polya-Hilbert conjecture which supposes 
that the imaginary part a of the zeros v = 1/2 + ia corresponds to eigenenergies of some 
unknown quantum mechanical hamiltonian. The well-studied statistical properties of the 
zeros, in particular the Montgomery- Odlyzko law and its relation to the gaussian unitary 
ensemble of random matrices 7, sj], led Berry to propose a quantum chaos interpretation of 
the oscillatory part of Riemann's counting formula for the zeros[9|. This is still within the 
context of quantum mechanics at zero temperature. Connes has proposed the complemen- 
tary picture that the zeros correspond to absorption lines|lo|. Sierra has recently proposed 
a consistent quantization of the Berry-Keating hamiltonian based on the Russian doll model 
of superconductivity, which possesses a cyclic renormalization group flow HI] . 

The fundamental duality relation relating ((u) to ((1 — v) can be understood as a con- 
sequence of a special modular transformation for the quantum statistical mechanics of free, 
relativistic, massless particles in v — d + 1 spacetime dimensions. This is explained in 
the appendix, since it is tangential to the main development of this paper. In section VII 
Riemann's ( function on the critical strip is related to the quantum statistical mechanics 
of non-relativistic, interacting fermionic gases in Id with a quasi-periodic 2-body potential 
which depends on v. This is thus a many-body problem at finite temperature. It is perhaps 
not completely unanticipated that the RH could have a resolution in the present context, 
since one of the very first places in physics where ( appeared was in Planck's work on black 
body radiation, which is widely acknowledged as the birth of quantum mechanics 12|. The 
black body theory is the bosonic quantum statistical mechanics of a 3c? gas of free relativistic 
photons 271] . In this case what appears is ((A) = 7r 4 /90. As we describe here, in order to 
get into the critical strip one needs d = 1 non-relativistic fermions interacting with a quasi- 
periodic potential. The fermionic nature renders the required integrals convergent and is 
also needed for a well-defined hamiltonian. This quasi-periodicity implies our approach is 
closest to Sierra's since there also a periodicity was important, but in the renormal- 
ization group at zero temperature. Our approach is thus essentially different just because 
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the physical context is different, however our interest in the RH actually stemmed from the 



work 



131 ] on cyclic renormalization group flows in Id relativistic systems and the observation 



that the finite temperature behavior was in part characterized by ((1 — ia), where 2ir/a 
is the period of the renormalization group flow[l^|. We now understand that one needs to 



consider non- relativistic models to get Re(v) < 1 in a natural way 28]. 

The understanding of the RH that emerges from the present context can be summarized 
as follows. When Re(v) > 1/2, there exists a two-body potential in position space that leads 
to a well-defined quasi-perodic kernel K in momentum space that is essentially a combination 
of Fourier transforms of the potential. The condition Re{y) > 1/2 comes about naturally 
as the condition for the convergence of the appropriate Fourier transform. The quantum 
statistical mechanics based on this kernel gives corrections to the pressure of the gas that are 
determined by a transcendental equation involving the polylogarithm Li u . To obtain this, 
one must work with the momentum-dependent kernel, i.e. the constant scattering length 
approximation made in earlier sections vanishes. If ((u) =0 then there would exist solutions 
with vanishing corrections to the pressure. Therefore the RH would follow from the simple 
physical property that non-zero interactions necessarily modify the pressure. We comment 
on the meaning of the actual zeros in section VII. 

II. FREE ENERGY AS A DYNAMICAL FUNCTIONAL OF FILLING FRAC- 
TION 

A. Generalities 

The free energy density (per volume) T is defined as 

T = --!- log Z, Z = Tt e -^ H -" N ) (2) 
P " 

where (3 = 1/T, p, are the inverse temperature and chemical potential, V is the d-dimensional 
spacial volume, and H and N are the hamiltonian and particle number operator. Since log Z 
is an extensive quantity, i.e. proportional to the volume, the pressure p of the gas is minus 
the free energy density: 

dhgZ 

V = T^ r = -T (3) 
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For most of this paper, we assume there is one species of bosonic (s = 1) or fermionic 
(s = — 1) particle. Given one can compute the thermally averaged number density n\ 

where k is the d-dimensional momentum. The dimensionless quantities / are called the 
filling fractions or occupation numbers. 

One can express T as a functional of / in a meaningful way with a Legendre transforma- 
tion. Define 

G = J-(fj) + fin (5) 

Treating / and fi as independent variables, then using eq. (j3J) one has that d^G = which 
implies it can be expressed only in terms of / and satisfies 5G/5f = fi. 

Inverting the above construction shows that there exists a functional F(f,fJ>) 

F(f,»)=G(f)-» J-0- d /(k) (6) 

which satisfies eq. flU) and is a stationary point with respect to /: 

The above stationary condition is to be viewed as determining / as a function of \i. The 
physical free energy is then T = F evaluated at the solution / to the above equation. We 
will refer to eq. ((7j) as the saddle point equation since it is suggestive of a saddle point 
approximation to a functional integral: 

Df e-WFW « e-PV* (8) 

In a free theory, the eigenstates of H are multi-particle Fock space states |k 1; k 2 ....). Let 
denote the one-particle energy as a function of momentum k. In this paper the theory 
is assumed to be non-relativistic with 

where m is the mass of the particle. It is well-known that the trace over the multi-particle 
Fock space gives 

™ = | / Wf log (1 " sc ~*" >) (10) 
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From the definition eq. (j3j) one finds the filling fractions: 

In order to find the functional F(f,fi) one first computes G from eq. (jHJ) and eliminates /i 
to express it in terms of / using eq. (jllj) . One finds 

Fo(/,A*) = / ^((wk-A*)/-^[(/ + a)log(l + a/)-/log/]) (12) 

One can then easily verify that SF /8f = has the solution / = / and plugging this back 
into eq. ( |T2l gives the correct result eq. ( ITUl) for JF . In the sequel, it will be convenient to 
trade the chemical potential variable \x for the variable fa: 

M/, ft) = 4 / ^ (• W + ./) + /log (^^) ) (13) 

There is another way to view the above construction which involves the entropy. Write 
eq. (TI2]) as 

F=S- l -S (14) 

where £ is the first (u — fi)f term in eq. (fl2l) . which is the energy density, and S is the 
remaining term in brackets. One can show by a standard counting argument, which involves 
the statistics of the particles, that S represents the entropy density of a gas of particles. 
(See for instance is|.) 

Let us now include interactions by writing 

F(/,/o) = Fo(/,/o)+Fi(/) (15) 

where Fo is given in eq. (1T3]) and we define U as the "potential" which depends on / and 
incorporates interactions: 

F ^-pjWf u{m) (16) 

Given F, f is determined by the saddle point equation SF /Sf = 0. It is convenient to define 
a pseudo-energy e as the following parameterization of /: 

(17) 

Then the saddle point equation and free energy density take the form: 

e = u-»---^ (18) 

T = -\ I Wf [ ~ s log(1 - se ~ P£) fdf)u] (19) 
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B. Two-body approximation 

It was shown in [l| how to express U in terms of S-matrix scattering amplitudes. In 
general there are terms involving the forward scattering of n — > n numbers of particles for 
n > 2. The most important is the two-body contribution. It takes the form 

tf(k) = |/(k) (K*/)(k) (20) 

where we have defined the convolution 

(K*/)(k) = |^K(k,k')/(k'), (21) 

and the kernel K is given by the 2-particle to 2-particle forward scattering amplitude: 

K(k,k / )=^i2;i2(k,k / ) (22) 

How to compute the kernel from the hamiltonian is described below. 

In terms of the pseudo-energy, the saddle point equation and free energy take the following 
forms: 

e = u-fi-K* ] (23) 



T = — : 



1 / d d k f . ^ if 



-s log(l - se~ pe ) + ^-(e - u) + //) 



(24) 



/3 J (2vr) d 

III. SECOND QUANTIZED HAMILTONIANS 

In this paper we consider second quantized hamiltonians of the form: 

H = J d d x (^W> f • VVM + ^Jd d *J V(x - x') ^ t (x)^ t (x / )^(x')^(x) (25) 

where V is the 2-body potential. The field satisfies the canonical commutation relations 

^(x)V + (x') - s ^(x'M*) = <J(x - x') (26) 

where again s = ±1 corresponds to bosons/fermions. Expanding the field in terms of 
annihilation operators, 

^ (x) = / eikx ak (27) 



S 



this leads to the canonical (anti-) commutation relations: 

a k al - s a£,a k = (27r) d 5(k - k') (28) 

The free Hilbert space is thus a bosonic or fermionic Fock space with momentum eigenstates 
normalized as: |ki, k 2 , k n ) = ■ ■ ■ aJ k |0). 

In this paper we only consider the lowest order contribution to the kernel K. Let H\ 
denote the interacting part of H which depends on the 2-body potential V. Then to lowest 
order, 

(a'lHtla) = -(2n) d 5 {d) (k a - k a >) M a >. a (29) 

where M. is the scattering amplitude for the asymptotic states \a). Therefore to lowest 
order the kernel is given by 

K(k,k') = -V- 1 (k,k'|#i|k,k') (30) 



IV. HARD CORE BOSON MODEL REVISITED 



In this section we consider bosonic particles with a delta-function two-body potential: 

V(x,x') =7<5 W (x-x') (31) 

The model has been solved exactly in Id by Lieb and Liniger UM , and the thermodynamic 
Bethe ansatz was first discovered in the context of this model[2|. 
Using (2vr) d (5( d )(0) = V, one finds 

K = - 7 (32) 

Positive 7 corresponds to repulsive interactions. Since the kernel is constant in this ap- 
proximation, it is equivalent to the constant scattering length approximation made in the 
literature; hence we do not obtain any new results for the hard core Bose gas here. We 
include this section and the next mainly for illustrative purposes and as a warm-up to the 
sequel. 

The coupling constant 7 has units of energy x volume. It can be expressed in terms of a 
physical scattering length a as follows. To first order in perturbation theory the differential 
cross-section in the center of mass is: 

(33) 



dQ 4(2tt) 



d-l 



where k is the magnitude of k for one of the incoming particles. Since a cross section has 
dimensions of length^" 1 , we define a such that the cross section is a d ~ x when the wavelength 
of the particle is 2n/a: 

% ~ " d " < 34 > 

ail k~2n/a 

This leads us to make the definition: 

7 a d ~ 2 

We carry out our analysis for arbitrary spacial dimension d > 0. Using J d d k = 
/ ^ ^ d_1 where T is the standard T-function, rotationally invariant integrals over 
momenta k can be traded for integrals over u^. 

«* '^^r^j^n (36) 



(27r) d V27T/ r(d/2) y o 
For a constant kernel K = —7, and d > 0, the solution to the integral equation eq. (f23j) 
takes the simple form: 

e(k) = cu k - + T6 (37) 
where 5 is independent of k. To determine 6 one needs the integral: 

dx = T(u)U v {z) for Re(u) > (38) 



The function \A v (z) is the standard polylogarithm, defined as the appropriate analytic con- 
tinuation of 

00 n 



n=l 



For fixed v, the function XA v {z) has a branch point at z = 1 with a cut along Re(z) > 1, 
and the above integral is not valid at z — 1. We will later need the fact that when z = 1, 
the above integral is only convergent when Re{v) > 1 and is given in terms of Riemann's ( 
function: 

roo v-l 

/ dx^— = r(u)C(v), Re{v)>\ (40) 
Jo e — 1 

Equivalent ly, 

Li„(l) = C(f), Re(u)>l (41) 
Eq. (1231) then leads to the following equation satisfied by 5: 

5 = h T Li d/2 (z^z s ) (42) 
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where we have defined the fugacities 



Zfl = e^, z 5 = e- 5 (43) 



and a renormalized thermal coupling hx and thermal wavelength A^: 



d-2 



To simplify subsequent expressions, we henceforth set the mass m — 1/2, unless otherwise 
stated. This leads us to define: 



T ^Af = £- (m=1 / 2) < 45) 



The equation ([4*2]) is a transcendental equation that determines 5 as a function of /j,, T, 
and the coupling h T . Given the solution S(fi) of this equation, using eq. (124]) the density 
can be expressed as 

n(fi) = f d ^U d/2 ( Zfl z 5 ) (46) 
For d > one can integrate by parts and obtain the following expressions for the free energy 



T = -T T d/2 ^Li (d+2 )/2(^<5) + ^Li d/2 (z M 2; 5 )^ 



(47) 



Bose-Einstein condensation can be described rather generally as follows. The property 
that an extensive, i.e. proportional to the volume, number of particles occupy the ground 
state with |k| =0 implies that the filling fraction / diverges at some critical chemical 
potential fi c : 

lim /(|k| = 0,/i) = oo (48) 



Given fi c , one can speak of a critical density n c : 

nc{T) = /(2^ /(k!/ic) (49) 
One also has a critical temperature T c defined as 

n c (T c ) = n phys , (50) 

where n p h ys . is the physical density. The reason that T c corresponds to a transition is that 
because of the divergence, one must treat the density of particles in the ground state n gs 
separately from the above expressions for n: 

n p h ys . = n gs + n (51) 
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where n gs is the density of particles in the ground state and n is the integral of /. When 
T > T c , n gs vanishes. The number of particles in the ground state doesn't actually vanish, 
but rather it is no longer proportional to the volume, so that n gs effectively vanishes. 
In the pseudo-energy description (|17|) of the filling fraction /, the condition (|48p is just 

e(|k| = 0,f/ = // c )=0 (52) 

Eq. (1231) then leads to the equation for /i c : \i c = T5(fi c ). In terms of the fugacities the 
critical point is the simple relation 

^c^(mc) = 1 ( 53 ) 
At the critical point, the arguments of the polylogarithms are 1, and if d > 2, then 1^/2(1) = 
((d/2). (Eq. dSD). This in turn implies that the critical temperature and density are 
independent of the coupling /t T : 

fjL e = h T ((d/2) (54) 
n c = ((d/2)T d / 2 (55) 

^ - 4 * fe) (56) 

Though the critical temperature and density do not depend on the coupling Ht to this order, 
the free energy at the critical point has corrections that do: 

^=-C(^)TT rf /^l + ^i^j (57) 

The ((u) function has a simple pole at v — 1, and this has a physical significance in 
the present context. In particular, it implies that the critical density n c is infinite and the 
critical temperature T c is zero in 2d to the order we have calculated. This means that in 
3d Bose particles have a stronger tendency to Bose-Einstein condense in comparison with 
2d. This simple pole at v — 1 can thus be interpreted as a manifestation of the Mermin- 
Wagner theorem, which states that finite temperature continuous phase transitions are not 



possible in 2d 18j . Another way of viewing this is that in 2d, bosonic particles behave more 
like fermions and don't readily Bose-Einstein condense. This fact is ultimately responsible 
for why one needs to treat bosonic particles with fermionic exclusion statistics in Id as far 
as their quantum statistical mechanics is concerned. See for instance 171 ] . For relativistic 
models also, it appears that only the fermionic thermodynamic Bethe ansatz equations are 
consistent |l9|, |2Cj . 
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V. FERMI GASES AND THE FERMI ENERGY 



In this section we consider fermionic particles with the same constant kernel K = —7 
as in the last section. To obtain such a kernel from the second quantized hamiltonian one 
needs at least two species of fermions since ip 2 — 0, however for simplicity we consider only 
one species. 

The integral we need is now: 

/ dx — = -T{v)U v {-z) Re{v) > (58) 

Jo e x + z 

It is important to realize that contrary to the bosonic case, eq. ()40l) . the above integral is 
valid when z = 1 throughout the critical strip < Re{y) < 1 and is given by the above 
formula where: 

- Li„(-1) = (1 - 2 1 - 1 0C(z/) Vu (59) 

In physical terms, fermions are much more stable than bosons in d < 2 because of the Pauli 
exclusion principle. 

The pseudo-energy still takes the form (1371) . and the saddle point equation (1231) now leads 
to the equation: 

5 = -h T \ji d / 2 {-z^z & ) (60) 
The density and free energy can be expressed as: 



n 



-T d / 2 U d/2 {-z^z s ) (61) 



T = T f d ' 2 {hi {d+2)/2 {-z^z 5 ) + -U d/2 (-z^z 5 ) 

In the fermionic case, the considerations of the last section are replaced by the concept 
of a Fermi surface. As we'll see, mathematically there are analogies between the existence 
of the Fermi surface and Bose- Einstein condensation. The Fermi energy ujf is defined as 
the uppermost energy that is occupied, i.e. has / 7^ 0. In the limit of zero temperature 
and no interactions, / is a step function: /(k) = 1 for |k| < |kp| and zero otherwise, where 
ujp = hp 2 /2m. At finite temperature the sharp step is broadened. 

The formulas eqs. fl6QII6ip can be used to study the effect of interactions and finite 
temperature on to p. The Fermi energy naturally can be defined as the point where / = 1/2, 
i.e. where e = 0: 

e(k F ) = = cup - fi + T5 (62) 
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The latter implies: 

Zli z s = e^ F (63) 
The Fermi energy may now be expressed in terms of the density n: 

n = -f d / 2 Li d/2 (-e^) (64) 

The above formula determines top as a function of T and n. As the temperature goes to 



zero, one can use 



lim U v {-e^' T ) « -fe^ (65) 

t^o vy ' r(z/+i) v ; 



to obtain 

2/d 



^ = 47r(r(^)n)' /a (66) 



In 3d the above formula is equivalent to well-known results 15]. 



Note that, to the order we have done the computation, the Fermi energy does not depend 
on the interactions since eq. AMI) no longer depends on the interaction. However in the 
presence of interactions 5^0, the free energy is modified from its non- interacting value, as 
eq. (IBTj) shows. 



VI. 2d BOSE AND FERMI GASES 

In this section we consider the 2d case which is rather special. Refering to eq. (144j) the 
coupling hx is temperature independent in 2c?, and is in fact dimensionless. This implies 
that in 2d the scattering length a d ~ 2 should simply be replaced by a dimensionless coupling 
we will denote as h. Many of the previous formulas are especially simple due to the identity 
hii(z) = — log(l — z). We also describe a strong similarity to the conformally invariant limit 
of relativistic systems in one lower dimension. 



A. Bosonic case 



The equation ff42|) for 8 becomes algebraic: 

z s = {l- z^f (67) 

The density is simply a logarithm 

T 

n = -— log(l - z^zs) (68) 
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and the free energy is expressed in terms of the dilogarithm: 

rp2 / r \ 

^ = ~4tt ( U ^ z ^ ~ 2 l0g ^ ~ Z ^ ) ^ 
The density is well defined and positive at zero chemical potential /i, so long as < z$ < 1. 
The equation fl67|) has a real solution in this range as long as h > 0, i.e. as long as the gas 
is repulsive. This observation means that the bosonic instability described in section IV 
can be cured by repulsive interactions. For the remainder of this section we set z^ — 1 and 
denote Zg simply as z. Let us define c as the following coefficient in the free energy: 

This is analagous to the relativistic case in Id, where one defines T = —c7iT 2 /Q where in 



that context c is the Virasoro central charge 23, |24|. In terms of the solution to the equation 



( I67p . c is given in terms of Roger's dilogarithm: 

c(z) = - 2 Lr 2 (z) (71) 

where 

Lr 2 (z) = U 2 (z) + \ log |z| log(l - z) (72) 

Interestingly, the above formula for c in terms of Lr 2 is identical to the 
when one studies the conformal limit of relativistic TBA systems in l<ij" 



The function Lr 2 is known to satisfy the following functional relations 



brmulas that arise 

B,y. 



Ti 2 



Lr 2 (^)+Lr 2 (l-^) = {(2) = — (73) 

o 

hi 2 {z) + Lr 2 (-z/(l - z)) = 
Using the first relation, one sees that Lr 2 (l/2) = £(2)/2. This, and the free fermion and free 



boson cases (h = 0) were known to Euler[2lJ. Landen found more relations as follows 22]. 
If r is a root to the polynomial equation z 2 + z = 1, then the above functional equations 
become linear equations for Lr 2 (z) with argument z = r, — r, — 1/r and r 2 . Let us chose 
r = (-\/5 — l)/2 which is the golden mean. When z is positive, the above relations imply 
that for special values of the coupling h, c is a rational number: 



(74) 



h = 


o. 


z 


= 1, 


c = 


1 


h = 


1/2, 


z 


= r, 


c = 


3/5 


h = 


1, 


z 


= 1/2, 


c = 


1/2 


h = 


2. 


z 


= r 2 , 


c = 


2/5 
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The above relations may be viewed as results in additive number theory by virtue of eq. 

Since the pressure p = —T, the coefficient c is a measure of the pressure of the gas. 
One observes from eq. (171]) that as h increases the pressure decreases. This is expected on 
physical grounds: larger h means stronger repulsive interactions so the gas is less dense. 

B. Fermionic case 

In this case the equation that determines z (at zero chemical potential) is 

z = (l + z)- h (75) 

The density takes the form 

n = £ log(l + z) (76) 
The free energy has the same form as in eq. (170]) where now 



o = ~ Lr 2 (-*) (77) 

Because of the Pauli exclusion principle for fermions, the fermionic gas is more stable 
than the bosonic one. First note that unlike the bosonic case, the eq. (175]) continues to 
have real solutions even if the interactions of the gas are attractive with h < 0. Another 
important feature is that whereas in the bosonic case Li2(z) has a branch cut along the real 
z axis from 1 to oo, there are no branch points for z along the negative axis. So for fermions, 
the free energy is well defined for any real z < 1. 

The special rational points of the bosonic case also have a fermionic version: 

h = — 1, z = oo, c = 1 

h = -1/2, z = 1/r, c = 3/5 (78) 

h = 0, z = 1, c= 1/2 

h = 1, z = r, c = 2/5 

Note here also that increasing h decreases the pressure since the gas is either less attractive 
or more repulsive. 
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C. Many mixed particles 



For far we have only considered a single bosonic or fermionic particle. It is straightforward 
to extend this to many types of particles of mixed statistics. Let m a and s a — ±1 denote 
the mass and statistics of the a-th particle. In the two-body approximation, we consider the 
following contribution to Fi' 

U = \ £>(k)(KW>)(k) (79) 

a,b 

where K. ab is the scattering amplitude of a with b types of particles. The saddle point 
equation now reads: 

e a = u a - n a ~y~] K ab * —~— (80) 

b 

For a constant kernel K. ab the above equation has the solution: 

£a = U a - fl a + T5 a (81) 

In 2d the equation satisfied by S a is again algebraic: 

Za = - s b z h ) h <*"> (82) 

b 

where z a = e~ Sa and h ab are dimensionless coupling constants that parameterize the kernel 
K a 5. The density n a of the a-th species is given by 

TTl T 

n a = -s a ^— log(l - s a z a ) (83) 
The total free energy is given by eq. ( 1701) with 

c = 2^2 m a c Sa {z a ) (84) 

a 

where c + is the bosonic expression (ITT]) and c_ the fermionic one (1771) . 

The exists a vast number of examples where certain choices of m a , s a , and h ab lead to 
rational c. For instance, one can translate known fermionic relativistic systems in Id to 
the present 2d non-relativistic context. Some of the latter are known to be related to root 

systems of Lie algebras^ (See for instance [3].) There are also many more examples that 

follow from results in [25|, not all of which have Id relativistic analog. For purposes of 
illustration, consider a two-particle theory with one boson and one fermion, of equal mass 
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1/2, and with h a b = 1, Va, b. This structure suggests a supersymmetric theory. The z\,z 2 
are solutions of 

21 = irt- 22 = ttS (85) 

which implies 24 = z 2 = — 1. It can be verified that this theory has c = 3/4. 



VII. QUASI-PERIODIC KERNELS AND THE RIEMANN HYPOTHESIS 

In order get into the critical strip of ((u) with < Re{y) < 1, inspection of our work so 
far suggests that one could try to analytically continue in the spacial dimension d to complex 
values. However it would remain unclear what complex d actually means physically. A more 
physical approach is to consider non-constant, quasi-periodic kernels in fixed dimension d 
an integer. As we will show, d — 1 is sufficient to cover the whole critical strip. Because 
of the previously discussed bosonic instabilities for d < 2, one must deal with a fermionic 
gas. In the next subsection we will simply hypothesize a certain kernel K and work out its 
consequences. In the subsequent subsection we will show how to obtain such a kernel from 
a second quantized hamiltonian with a 2-body potential that is also quasi-periodic. 



A. Quasi-periodic kernel in Id 

We assume that the 2-body potential V(x, x') in eq. ( 1251) is translationally invariant, i.e. 
it depends only on the difference x — x'. In momentum space this implies that the kernel 
also depends only on the difference: 

K(k,k') =K(k-k') (86) 

We also assume rotational invariance so that K depends only on k = |k|. Let us suppose 
there exists a hamiltonian which leads to the following kernel: 

K(k) = -Re^k 2 "- 1 ) (87) 

where v is assumed to be a complex number and ^ v is a constant. Note that the lowest 
order constant scattering length approximation vanishes as long as Reiv) > 1/2. 
The saddle point equation, with zero chemical potential, reads: 

r dv i 
<*) = ^ + J^ Re W k - k 'l 2 ^) ^oTT (88) 
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At low temperature, e(k) can be approximated near k = 0, where it continues to take the 
form eq. fl37|) . with fi = 0. The constant 5 must satisfy the equation: 

5 = -Re \T u ~ l h u U v (-z s )] (89) 

where h u = 7^r(z/)/27r and as before zs = e~ s . The density and free energy however have 
the same expressions as in section IV specialized to d — 1: 

n = -{T/A-k) 1 ' 2 Li 1/2 (-z«) (90) 
T = (T 3 /47r) 1/2 (u 3/2 (-z s ) + 5 -Li 1/2 (-z s )\ 

There are some trivial 5 = solutions to eq. ( |89l) which arise when the pref actor 1 — 2 l ~ u 
in eq. (1591) equals zero. We can remove them by chosing 7^ = ((1 — 2 1 ~ u )r(v))~ 1 . We have 
divided by T(z/) in order to remove another trivial zero which can arise from r(zoo) = 0. 
The coupling constant then becomes: 

o("-3)/2 

hu = 7 r (9!) 

2 7 rsinh^(i=^ 2 ) X 



2 

We can now give a clear meaning to a zero v of the zeta function ((u) = 0. Because of 
the relation fl5T?|) . when v is a zero, then a solution to eq. (I8T?|) is 5 = 0. This solution exists 
for an?/ temperature T. By plotting the left and right hand sides of eqn. (I8"9"|) . one sees that 
in general there is another solution at 5 7^ 0, but the latter depends on T, and in fact for 
some T this solution disappears. Thus for the generic 5 = solution, this means there are 
no corrections to the free energy at arbitrary temperature T, i.e. the density and free energy 
are the same as in a free (non-interacting) theory: 

n = (T/4tt) 1/2 (1-V2)C(1/2) (92) 
T = -(T 3 /4vr) 1/2 (1-1/^(3/2) 

Note that the density is still positive since C(l/2) is neg ative: (1 - V2)((l/2) w .60649. The 
pressure is also positive: (1 - l/ v / 2)C(3/2) w .76515. 

The Riemann hypothesis would then follow from the line of reasoning: If (i) The leading 
contributions to the pressure of the gas are the ones calculated in this work, (ii) There is 
a hamiltonian that leads to the //-dependent quasi-periodic kernel (187j) . and (iii) Non-zero 
interactions necessarily modify the pressure of the gas over a range of temperatures, then 
can have no zeros. In the next subsection we address (ii). 
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B. Real space potentials 



We now show that there are indeed real space hamiltonians in Id that lead to the above 
quasi-periodic kernel. Let v (k) be the Fourier transform of the 2-body potential: 

V(x-x') = J -0- d e^—') v(k) (93) 

Using eq. ( 1301) one finds 

K(k) = --(u(k)+u(-k) + 2su(0)) (94) 

We have included the statistics parameter s in order to point out some features. 
Consider the following potential: 

= Re (^) (95) 



n |x| 

where b v is a constant. If b v is real and v = a/2 + iot then the potential is quasi-periodic in 
log \x\: 

V(x) = 6 - C ° S y x|) (96) 

This model is not known to be integrable, unlike the case of a delta-function potential, thus 
there is no known exact thermodynamic Bethe ansatz. The lowest order corrections are the 
ones calculated in the last section. 

When the particles are fermionic with s = — 1 the v (0) term in eq. (|9~4"|) precisely renders 
the kernel well behaved as k — > 0. The kernel would be singular if the particles were bosonic. 
The kernel is then given by the following integral: 

K(k) = Re (2b v J dx x~ 2u sin 2 (A;x/2)^ (97) 

This integral is convergent if the following condition is satisfied: 

1/2 < Re(u) < 3/2 (98) 

If the above condition is met, the kernel is the following: 

K(Jfe) = -Re [b u k 2v ~ x sin(7rz/)r(l - 2v)\ (99) 

Finally we can chose b v in order to obtain the kernel in eq. (19T|) . Using the identity 

sin(7rz/)r(l - 2v)Y{v) = ^2^(1/2 - v) (100) 
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this fixes b v to be: 



b, 



2(5i/-3)/2 



(101) 



v 



0Fr sinh ( ci-^iogeo 



) 



The condition (!98|) is precisely what one needs for the RH. What is then the meaning of 
the known zeros at Re{y) = 1/2? These are models with the kernel fl87|) that give vanishing 
leading corrections to the pressure. The calculations of this section show that such a kernel 
does not arise in a convergent manner from a real space potential since the condition (|98p 
is violated. Note that the T function in eq. (1991 develops a pole at v — 1/2, which suggests 
that an additional low-energy regularization could still lead to sensible models with the 
kernel flH7|) that provide physical realizations of the Riemann zeros. 

VIII. CONCLUSIONS 

We have shown how the formalism developed in [l] can lead to new results for the quantum 
statistical mechanics of interacting gases of bosons and fermions. Our main results were 
summarized in the introduction. 

Clearly our most interesting result is the formulation of the Riemann hypothesis in the 
present context. We have essentially given a physical argument from which it follows. In 
order to develop this argument into a rigorous mathematical proof, one mainly needs to give 
a firm foundation to the theoretical methods used here, namely the framework developed 
in [l|. One also needs to understand more rigorously how the contributions considered in 
this paper are the leading ones, i.e. one needs clearer control of the approximations we have 
made. 
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X. APPENDIX 

In this appendix we explain how the v — > 1 — v duality of Riemann's zeta function ((u) 
can be understood as a special modular transformation in a Lorentz-invariant theory. 
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Consider a free quantum field theory of massless bosonic particles in d + 1 spacetime 
dimenions with euclidean action S = J d d+1 x (<90) 2 . The geometry of euclidean spacetime 
is take to be S 1 x R d where the circumference of the circle S l is (3. We will refer to the S 1 
direction at u x" . Endow the flat space R d with a large but finite volume as follows. Let us 
refer to one the directions perpendicular to x as y with length L and let the remaining d—1 
directions have volume A. 

Let us first view the x direction as compactified euclidean time, so that we are dealing 
with finite temperature T — 1/(5. As a quantum statistical mechanical system, the partition 
function in the limit L, A — * oo is 

Z = e-w™ (102) 
where V = L ■ A and T is the free energy density. Standard results give: 

r( ' 3) = J) I Wf log (1 " e ~ m) (103) 

The euclidean rotational symmetry allows one to view the above system with time now 
along the y direction. In Id, interchanging the role of x and y is a special case of a modular 
transformation of the torus. In this version, the problem is a zero temperature quantum 
mechanical system with a finite size (3 in one direction, and the total volume of the system 
is V — P • A. The quantum mechanical path integral leads to 

Z = e - LEo ^ (104) 

where E is the ground state energy. Let £ = E /V' denote the ground state energy per 
volume. Comparing the two "channels", their equivalence requires S ((3) = J-{(3)- In this 
finite-size channel, the modes of the field in the x direction are quantized with wave-vector 
k x = 2irn/{3, and the calculation of £ is as in the Casimir effect: 

f » = ^Z/(^r( k2 + ( 2 ™/« 2 ) ,/2 (105) 

The free energy density T can be calculated using / d d k = 2-K d / 2 /T{d/2). For d > the 
integral is convergent and one finds 

i r(d+i)c(d+i) . , 

p d +i 2 d -H d ' 2 Y(d/2)d K } 
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For the Casimir energy, So involves ^ngg \ n \ d which must be regularized. As is usually 
done, let us regularize this as 2£(— d). Then: 

So = ^ 2 r(-rf/2)C(-rf) (107) 

Define the function 

£M = ^ 2 T(u/2)C(u) (108) 
Then the equality £ = T requires the identity: 

£(1-1/) (109) 

The above relation is a known functional identity that can be proven using complex analysis. 
Thus we have demonstrated that ( function regularization of the Casimir energy is consistent 
with a modular transformation to the finite-temperature channel. On the other hand, our 
calculations can be viewed as a proof of the identity (jlU9p based on physical consistency. 



[1] A. LeClair, Quantum statistical mechanics of gases in terms of dynamical filling fractions and 



[2] 
[3] 
[4 
[5 
[6] 
[7] 
[8] 
[9] 

[10 

[11 
[12 

[13 



scattering amplitudes, J. Phys. A 40 (2007) 9655 arXiv:hep-th/0611187 



C. N. Yang and C. P. Yang, Jour. Math. Phys. 10, (1969) 1115. 

R. Dashen, S.-K. Ma and H. J. Bernstein, Phys. Rev. 187 (1969) 345. 

T. D. Lee and C. N. Yang, Phys. Rev. 113 (1959) 1165; Phys. Rev. 117 (1960) 22. 

Al. Zamolodchikov, Nucl. Phys. B342 (1990) 695. 

H. M. Edwards, Riemann's Zeta Function, Academic Press, New York (1974) 
H. Montgomery, Analytic Number Theory, AMS (1973). 
A. Odlyzko, Math. Comp. 48 (1987) 273. 

M. V. Berry, in Quantum Chaos and Statistical Nuclear Physics, eds. T. H. Seligman and H. 
Nishioka, Lecture Notes in Physics, No. 263, Springer Verlag, New York (1986). 
A. Connes, Selecta Mathematica, New Series 5 (1999) 29, j arXiv:math.NT/9811068| 



G. Sierra, J. Stat.Mech. 0512 (2005) P006, arXiv:math.NT/0510572 



M. Planck, Annalen der Physik, vol. 4 (1901) 553. 



A. LeClair, J. M. Roman, and G. Sierra, Nucl.Phys. B675 (2003) 584, hep-th/0301042; A. 



LeClair and G. Sierra, J.Stat.Mech. 0408 (2004) P004, hep-th/0403178. 

23 



A. LeClair, J. M. Roman, and G. Sierra, Nucl.Phys. B700 (2004) 407, |hep-th/0312141 



L. D. Landau and E. M. Lifshitz, Statistical Physics, Pergamon Press (1980). 
E. Lieb and W. Liniger, Phys. Rev. 130 (1963) 1605. 



V. Gurarie, Phys. Rev. A73 (2006) 033612, cond-mat/05 11411 



N. D. Mermin and H. Wagner, Phys. Rev. Lett. 17 (1966) 1133; P. C. Hohenberg, Phys. Rev. 
158 (1967) 383. 

T. Klassen and E. Melzer, Nucl. Phys. B338 (1990) 695. 



G. Mussardo and P. Simon, Nucl. Phys. B578 (2000) 527, |hep-th /9903072, 

L. Euler, Institutiones calculi integralis, 1 (1768) 110-113. 

J. Landen, Mathematical memoirs, 1 (1780), 112. 

J. Cardy, J. Phys. A17 (1984) 385. 

I. Affleck, Phys. Rev. Lett. 56 (1986) 746. 

L. Lewin, ed. Structural Properties of Polylogarithms, Mathematical Surveys and Monographs, 
vol. 37, 1991. 



For a comprehensive list see M. Watkins at http:/ /secamlocal.ex.ac.uk~mwatkins/zeta/physics.htm 



It is not clear from Planck's paper whether he was aware he was dealing with Riemann's £ 
function. He simply writes l + p + ^ + p + -- - = 1.0823 
[28] There are also some convergence problems with the usual TBA, which are avoided in the 
present work. 



24 



